Finite-temperature phase transition in a class of 4-state Potts antiferromagnets 
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We argue that the 4-state Potts antiferromagnet has a finite-temperature phase transition on any 
Eulerian plane triangulation in which one sublattice consists of vertices of degree 4. We furthermore 
predict the universality class of this transition. We then present transfer-matrix and Monte Carlo 
data confirming these predictions for the cases of the union-jack and bisected hexagonal lattices. 
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The (/-state Potts model 0, Q plays an important 
role in the theory of critical phenomena, especially in 
two dimensions [iMa], and has applications to various 
condensed- matter systems . Ferromagnetic Potts mod- 
els are by now fairly well understood, thanks to universal- 
ity; but the behavior of antifcrromagnctic Potts models 
depends strongly on the microscopic lattice structure, so 
that many basic questions must be investigated case-by- 
case: Is there a phase transition at finite temperature, 
and if so, of what order? What is the nature of the low- 
temperature phasc(s)? If there is a critical point, what 
are the critical exponents and the universality classes? 
Can these exponents be understood (for two-dimensional 
models) in terms of conformal field theory [H[? 

One expects that for each lattice C there exists a value 
q c {£) [possibly nonintcger] such that for q > q c (C) the 
model has exponential decay of correlations at all tem- 
peratures including zero, while for q = q c {£) the model 
has a zero-temperature critical point. The first task, for 
any lattice, is thus to determine q c . 

Some two-dimensional (2D) antiferromagnetic mod- 
els at zero temperature can be mapped exactly onto a 
"height" model (in general vector- valued) @, Q- Since 
the height model must either be in a "smooth" (or- 
dered) or "rough" (massless) phase, the corresponding 
zero-temperature spin model must either be ordered or 
critical, never disordered. Experience tells us that the 
most common case is criticality The long-distance 
behavior is then that of a massless Gaussian with some 
(a priori unknown) "stiffness matrix" K > 0. The crit- 
ical operators can be identified via the height mapping, 
and the corresponding critical exponents can be pre- 
dicted in terms of K. Height representations thus provide 
a means for recovering a sort of universality for some (but 
not all) antifcrromagnctic models and for understanding 
their critical behavior in terms of conformal field theory. 

In particular, when the g-state zero-temperature Potts 
antiferromagnet on a 2D lattice C admits a height rep- 



resentation, one ordinarily expects that q = q c (C). 
This prediction is confirmed in most heretofore-studied 
cases: 3-state square-lattice JE S EH G3 , 3-state kagome 



13l |14| . 4-state triangular 1151, a nd 4-state on the line 
graph of the square lattice [141 Il6l | . The only known ex- 
ceptions are the triangular Ising antiferromagnet 12] and 
the 3-state model on the diced lattice (loj . 

Moore and Newman [l5[ observed that the height map- 
ping employed for the 4-state Potts antiferromagnet on 
the triangular lattice carries over unchanged to any Eule- 
rian plane triangulation (a graph is called Eulerian if all 
vertices have even degree; it is called a triangulation if all 
faces are triangles). One therefore expects naively that 
q c = 4 for every (periodic) Eulerian plane triangulation. 

Here we will present analytic arguments suggesting 
that this naive prediction is false for an infinite class of 
Eulerian plane triangulations, namely those in which one 
sublattice consists entirely of vertices of degree 4. More 
precisely, we predict that on these lattices the 4-state 
Potts antiferromagnet has a phase transition at finite 
temperature (so that q c > 4); we shall also predict the 
universality class of this transition. We will conclude by 
presenting transfer-matrix and Monte Carlo data con- 
firming these predictions for the cases of the union-jack 
[£>(4,8 2 )] and bisected hexagonal [D(4,6,12)] lattices. 

Exact identities. Let G = (V,E) and G* = (V*,E*) 
be a dual pair of connected graphs embedded in the plane 
(Fig. Hk). Then define G = {V U V* , E) to be the graph 
with vertex set V UV* and edges ij whenever i £ V lies 
on the boundary of the face of G that contains j £ V* 
(Fig. [TJa) . The graph G is a plane quadrangulation: on 
each face of G, one pair of diametrically opposite vertices 
corresponds to an edge e £ E and the other pair corre- 
sponds to the dual edge e* £ E*. In fact, G is nothing 
other than the dual of the medial graph M{G) = M{G*) 
[19^ . Conversely, every plane quadrangulation G arises 
via this construction from sonic pair G . 

Now let G be the graph obtained from G by adjoin- 
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Figure 1: (a) A dual pair G (black dots and lines) and G* 
(gray dots and lines), (b) The quadrangulation G (black/gray 
dots and thick lines) and triangulation G (all dots and lines). 

ing a new vertex in each face of G and four new edges 
connecting this new vertex to the four corners of the face 
(Fig. [TJd). This graph G is an Eulerian plane triangu- 
lation, with vertex tripartition V U V* U C where C is 
the set consisting of the "new" degree-4 vertices. Con- 
versely, every Eulerian plane triangulation in which one 
sublattice consists of degree-4 vertices arises in this way. 

We will show elsewhere [2(| that the 4-state Potts an- 
tiferromagnet at zero temperature (= 4-coloring model) 
on G can be mapped exactly onto the 9-state Potts 
/erromagnet on G and G*\ 

Z~(4,-l) =4-3-l y lz G (9,3) = 4-3-l y *lz G .(9,3) (1) 

where Zc(q,v) denotes the Potts-model partition func- 
tion with v — e J — 1 (J = nearest-neighbor coupling). 
The proof passes either via an RSOS model on G or via 
a completely packed loop model on A4(G). 

Height representation fldi l. Consider the 4-coloring 
model on an Eulerian plane triangulation 0. We can 
orient the edges of O such that the three edges around 
each face define a cycle (clockwise on one sublattice of O* 
and counterclockwise on the other). Let eo, ei, e.^ be unit 
vectors at angles 0, 2tt/3, 4tt/3 in the plane. Then, to any 
proper 4-coloring a of the vertices of 0, we assign heights 
h, in the triangular lattice such that h, — h; = eo, ei, e2 
on an oriented edge ij according as {o~i,o~j} = {1,2} or 
{3,4}, {1,3} or {2,4}, {1,4} or {2,3}. 

Phase transition and universality class. Let now G, 
G* and G be infinite regular lattices. Conformal field the- 
ory [H tells us that a g-state Potts ferromagnet with q > 4 
cannot have a critical point. Therefore the 9-state Potts 
ferromagnet in ([1]) is noncritical, suggesting that the 4- 
state Potts antiferromagnet at zero temperature on G is 
also noncritical [HI]. But since this model has a height 
representation, it cannot be disordered; therefore it must 
be ordered. It follows that the 4-state Potts antiferro- 
magnet on G has an order-disorder transition (whether 
first-order or second-order) at finite temperature. 

We can also understand the type of order in the 4- 
coloring model on G, and hence the universality class of 
the order-disorder transition in case it is second-order. If 
the lattice G is self-dual, the point (q, v) = (9, 3) lies on 
the self-dual curve v = yfq, which is expected to be the lo- 
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Figure 2: (a) Union-jack lattice, L — 6. (b) Bisected hexag- 
onal lattice, L — 8. The shaded areas show the minimal unit 
cells (pink) and the rectangular unit cells used in the row-to- 
row transfer-matrix computations (blue). The tripartition of 
the vertex set is shown in black/gray/white as in Fig. [T] 

cus of first-order transitions; so there are phases in which 
G is ordered and G* is disordered, and vice versa (nine 
of each) . We therefore predict that the 4-coloring model 
on G has phases in which the sublattice V is ordered in 
one of the four possible directions while V* and G are 
disordered, and the same with V and V* interchanged. 
The symmetry is S 4 x Z 2 , so we expect that the transi- 
tion is in the universality class of a 4-state Potts model 
plus an Ising model (decoupled). On the other hand, if 
G is not self-dual, then we expect (barring a fluke) that 
(q, v) = (9, 3) does not lie on a phase-transition curve; 
hence one of the lattices G, G* will be ordered (say, G) 
while the other is disordered. In this case we predict that 
the 4-coloring model on G has phases in which the sub- 
lattice V is ordered in one of the four possible directions 
while V* and G are disordered. The symmetry is S4, and 
the universality class is that of a 4-state Potts model. 

We recall that the central charge c and magnetic and 
thermal exponents X m ,X t are (c, X m , X t ) = (|, |, 1) for 
the Ising model and (1; g, 5) for the 4-state Potts model. 

Union-jack (UJ) lattice. The simplest self-dual case is 
G = G* = G = square lattice and G = union-jack lattice 
(Fig. [2^i). We computed transfer matrices with fully 
periodic boundary conditions for even widths L < 20 
(v = —1) [22| and L < 16 (general v). Estimates of 
c,X m ,X t were extracted from the free energy /l and 
free-energy gaps A/j, via 0] 

h = foe - W(6£ 2 ) + o(l/L 2 ) (2) 
Af L = 2ttX/L 2 + o(l/L 2 ) (3) 

Fig. [3] (upper left) shows estimates of c(v) at q = 4. 
The maximum at v ~ —0.95 indicates the transition: 
finite-size scaling (FSS) yields v c = —0.944(5) and c = 
1.510(5), in agreement with our prediction c = I + 3 = |- 
The crossings ofX m (v) and X t (v) yield v c = -0.9488(3), 
X m = 0.1255(6) and Xt = 0.51(2), in agreement with 
X m = I and X t = \ 

A similar plot for c(q) at v = — 1 shows the lattice- 
independent Bcrkcr-Kadanoff phase [c = 1 — 6(t— I) 2 /t 
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Figure 3: Estimates for central charge c and critical expo- 
nents X m , Xt for the U J lattice, as a function of v at q — 4 
(left) and as a function of q at v = —1 (right). Dashed vertical 
lines indicate our best FSS estimates of v c , qo and q c . Fits 
used Eqs. ©/© with o(l/L 2 ) replaced by A/L 4 , for three 
(resp. two) consecutive values of L. 



with q = 4cos 2 (7r/i)] for < q < qo and a noncritical 
phase for go < q < qc- The maxima of c(q) [Fig. [31 upper 
right] yield the estimates qo = 3.63(2), c = 1.43(1) and 
q c = 4.330(5), c = 1.63(1). The crossings of X m (q) and 
X t (q) yield q Q = 3.616(6), X m = 0.0751(3), X t = 0.88(2) 
and q c = 4.326(5), X m = 0.134(2), X t = 0.52(3). 

The data for go are consistent with qo = Bio = (5 + 
V5)/2 ps 3.61803 Q and c = 7/5, X m = 3/40, X t = 
7/8; the underlying conformal field theory could be a 
pair of m = 4 minimal models [25[ • 

Concerning q c , we have seen that at (q, v) = (4, v c ) the 
critical behavior is that of a 4-state Potts model plus an 
Ising model (decoupled), and it is compelling to think 
that this behavior will persist along a curve in the (q, v)- 
plane passing through (q,v) = (q c , — 1). However, it is 
possible that (q c , — 1) might be the endpoint of this curve, 
in which case the model could be driven there to some 
sort of multicritical behavior: for instance, a 4-state Potts 
model plus a tricritical Ising model (decoupled), which 



would have c : 
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0.13125 [25j. Alternatively, the critical exponents along 
the transition curve may vary continuously with q. 

We also simulated the q = 4 model, using a cluster 
Monte Carlo (MC) algorithm, on periodic L x L lattices 
with 8 < L < 512. We measured the sublattice mag- 
netizations Ma, -Mb, -Mc the nearest-neighbor correla- 
tions £ab,£ac,£bc and the next-nearest-neighbor cor- 
relations £aa,£bb,£cc- We then computed the 3x3 
sublattice susceptibility matrix and the 6x6 specific- 
heat matrix; from their eigenvalues we can extract the 
magnetic and thermal critical exponents. The leading 
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Figure 4: Monte Carlo data for the UJ lattice at q — 4. 
(a) Leading susceptibility eigenvalue Ai(x) divided by 
I/ 7//4 . (b) Leading specific-heat eigenvalue \i(C) divided 
by L. (c) Second susceptibility eigenvalue A2(x) divided by 
L 3 / 2 (log L)~ 1,/4 . Lines are meant only to guide the eye. 



susceptibility eigenvalue diverges with the predicted ex- 
ponent ~f/i> = 2 - 2X m = 7/4 (Fig. Hi), and FSS anal- 
ysis yields the estimate v c = —0.9485(1). Likewise, the 
leading specific-heat eigenvalue diverges with exponent 
a/v = 2 - 2X t = 1 [i| (Fig. gj>), and FSS analysis 
yields the estimate v c = —0.9483(2). It is curious that 



we do not see here the multiplicative logarithms that 
are expected [26| for the 4-state Potts model. The sec- 
ond susceptibility eigenvalue diverges as L 3 / 2 , probably 
with a multiplicative logarithm (Fig. 3J:), while the sec- 
ond specific-heat eigenvalue tends to a finite constant; we 
have no theoretical understanding of these behaviors. 

A transition in this model was recently predicted 
by Chen et al. [13, who found v c = -0.9477(5) by a 
tensor renormalization-group method; they also gave an 
entropy-counting argument predicting the type of order. 
However, in their approximation the specific heat was 
non-divergent, exhibiting a jump discontinuity. 
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Bisected hexagonal (BH) lattice. The simplest non- 
sclf-dual case is G = triangular lattice and G* = hexag- 
onal lattice, yielding G = diced lattice and G = bisected 
hexagonal lattice (Fig.^Ja). We computed transfer matri- 
ces with fully periodic boundary conditions for the same 
widths as for the UJ lattice, except that L must now be a 
multiple of 4 to be compatible with the periodic bound- 
ary conditions (see Fig. [2)3). FSS analysis of c(v) at q = 4 
yields the estimates v c = -0.8281(1) and c = 1.000(5), 
in agreement with our prediction c = 1. The crossing of 
X m (v) and the minimum oi X t (v) yield v c = —0.8280(1), 
X m = 0.15(1) and X t = 0.65(10), which are compatible 
with X m = i and A t = i although rather imprecise. 

The maximum of c(g) yields q c = 5.395(10), c = 
1.20(5). The crossing of X m (q) and the minimum of 
X t (q) yield q c = 5.397(5), X m = 0.15(1), X t = 0.6(1). 

We also did MC simulations for q = 4 and 8 < L < 512. 
The leading susceptibility eigenvalue diverges as expected 
as L 7 / 4 [possibly with a multiplicative (logL) -1 / 8 ] and 
yields v c = —0.828066(4). The leading specific-heat 
eigenvalue is compatible with the 4-state Potts behav- 
ior L(logL)~ 3 / 2 . 
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A transition in this model was also conjectured in [27 1. 

Our result q c > 5 suggests that there will be a finite- 
temperature transition also in the 5-state model. Quite 
surprisingly, we find this transition to be second-order, 
despite the absence of an obvious universality class (since 
q > 4); however, it is also conceivable that the tran- 
sition is weakly first-order, with a correlation length 
that is finite but very large. Preliminary results from 
transfer matrices are v' c = —0.9513(1), c = 1.17(5), 
X m = 0.16(1), X t = 0.56(4). Preliminary MC results 
are v' c = -0.95132(2), X m = 0.113(4), X t = 0.495(5). 
More detailed data will be reported separately 28 1. 

Taking into account the likely corrections to scaling, 
our data for (g, v) = (4, v c ), (5, v' c ) and (q c , —1) are com- 
patible with all three models being in the 4-state Potts 
universality class. 

Conclusion. We have given: (a) an analytical existence 
argument for a finite-temperature phase transition in a 
class of 4-state Potts antiferromagnets; (b) a prediction 
of the universality class; (c) large-scale numerics, using 
two complementary techniques, to determine critical ex- 
ponents; (d) determination of go and g c as well as v c ; 
and (e) the surprising prediction of a finite-temperature 
phase transition also for q = 5 on the BH lattice [28j |. 
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